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ABSTRACT 
We present mid-infrared spectra of six FeLoBAL QSOs at 1 < z < 1.8, taken with the Spitzer 
space telescope. The spectra span a range of shapes, from hot dust dominated AGN with sihcate 
emission at 9.7/im, to moderately obscured starbursts with strong Polycyclic Aromatic Hydrocarbon 
(PAH) emission. The spectrum of one object, SDSS 1214-0001, shows the most prominent PAHs yet 
seen in any QSO at any redshift, implying that the starburst dominates the mid-IR emission with 
an associated star formation rate of order 2700 M0 yr~^. With the caveats that our sample is small 
and not robustly selected, we combine our mid-IR spectral diagnostics with previous observations to 
propose that FeLoBAL QSOs are at least largely comprised of systems in which (a) a merger driven 
starburst is ending, (b) a luminous AGN is in the last stages of burning through its surrounding dust, 
and (c) which we may be viewing over a restricted line of sight range. 
Subject headings: galaxies: active - quasars: absorption lines - infrared: galaxies - galaxies: evolution 
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1. INTRODUCTION 

It is now well-established that there exist intimate 
links between the growth of supermassive black holes 
and the production of stars in galaxies. Indirect evi- 
dence for such links comes from, for example, the tight 
relationship bet ween central black hole and stellar bulge 
mass es fe.g. IMagorrian et al.l |T998l : IGebhardt et al.l 
120001 IGadotti fc Kauffmannl 120091) . and the coeval stel- 
lar populations in many massive elliptical galaxie s 



punlop et al.|[T99l lEUis et al.|[l997l: iRakos et al.i[2007l ) 



suggesting their stars formed within less than a Gyr 
of each other. Direct evidence comes from the ex- 
istence of galaxies that simultaneously harbor both 
high rates of star formation (of order tens to thou- 
sands Mf;) yr~^) and rapid accretion onto a central 
black hole (iSand ers fc Mirabell[l996l: iGenzel et al.l ll998l: 
iFarrah et alT l2003: LonsdaleeLalJ 120060 , all enveloped 
in large quantities of dust, making them extremely lu- 
minous in the infrared. The importance of these links 
for understanding the assembly history of galaxies over 
at least a substantial fraction of the age of the Uni- 
verse is demonstrated by the strong evolution of the lu- 
minosity function of IR-luminous galaxies with redshift 
(jSaunders et al.lll990HLe Floc'h et al.ll2005[ ). and the ex- 
istence of a profusion of IR-lumin ous sources in th e high- 
redsh ift Universe (e.g. iBarger et al. 1998: Eales et al.l 
ll999HCoppin et al.ll2006l: lAustermann et al.ll2009D . 
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An insightful way to study these links is to identify ac- 
tive galaxies at 'key' points, where they are either rapidly 
building up stellar or central black hole mass, or shifting 
from one evolutionary phase to the next. One such point 
may be the "youthful" QSOs; if the starburst precedes 
or is coeval with the QSCo, then young QSOs should be 
those in which the starburst is drawing to a close, and 
a QSO is starting to emerge from its surrounding dust. 
Such objects would provide an excellent laboratory for 
testing feedback mechanisms between the starburst and 
AGN (e.g. Silk & Rccs 1998; Ciotti & Ostriker 200;S 



Lagos et al.l 120081: iMoe et al.l 120091: iCeverino fc Klvvui 
20m . 



Finding the youthful QSO pop ulation is a challeng- 
ing and ongoing problem (e.g. ISanders et al.l 119881: 



Coppin et al.l 



Canalizo & Stockt o^ I200TI: iLacvl 120061: ------- 

20081 IGeorgakakis et al.l 120091 : iLipari et al.l l2009t )" 

nificant attention has fo cuse d on the Broad 

tion L ine (BAL) QSOs JLy nds 19 671: iWevmann et alj 
19911: iBrotherton et all [1998; .Schmidt fc HinesI |199£ 



Sig- 
Absorp- 



Arav et al.l 2001 



._ _ iGreenet al.l 12 001: Hah et al. 200L 

Reic hard et al.l l2003t iPriddev et al.i i2007t iGibson et ail 
2009), whose propertied^"! have been explained either 
as arising from being observed at a particular orien- 
tation, or because they are young objects still par- 
tially surrounded by dust. The properties of HiBALs 
are now thought to be primaril y an orientation ef- 
fect (|Surdei fc HutsemekersI 119871: iMurrav et all 119951 : 



^ The (optical) QSO phase is not necessarily the period in which 
the SMBH gains the bulk of its mass; indeed, there is evidence that 
the major periods of BH growth are coeva l with or precede the 
starbu r st, see e.g. IMartinez-Sansigre et al.l 120051) : IMerloni et al.l 
ll2010l) : nVeIster et al.l ll2oTof) 

^^ BAL QSOs come in 3 subtypes. High Ionization BAL QSOs 
(HiBALs) show absorption in CIV A1549, NV A1240, SilV A1394 
and Lya, and comprise about 10% of all QSOs (Trump et al. 2006). 
Low Ionization BAL QSOs (LoBALs) additionally show absorption 
in Mgll A2799 and other lower ionization species, and comprise 
~ 1.3% of all QSOs. Finally, FeLoBAL QSOs, in addition to show- 
ing all the absorpti on lines s een in LoBALs, also show weak iron 
absorption features IHazard' et al. 1987; Becker et al. 1997). They 
account for around 0.33% of all QSOs, though the exact fraction 
is unclear. 
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Schmidt fc HhleifTQQl [Gallagher et aP [20071: [DoreTall 
20091) ■ while the LoBALs remain controversial, with 



evidence f avoring both orientation and evolutio n (e.g 



Voit et all 1993; Ogle ct al. 1999; GaUagher et all 



2007 : 



Ghosh fc P unslv 2007; Montenegro- Montes et all 12008 : 
Urrutia et al. 2009 ; Zhang et al. 2010l ). For FeLoBALs 



there is also controversy, but the evidence favoring the 
evolution scenario is (arguably) stronger than for the 
LoBALs. Based on rest-frame UV and optical spec- 
tra. Hall et al (2002) suggest that FeLoBALs are young 
objects still enveloped in dust. Similar conclusions are 
reached by Gregg et al (2002), who also postulate that 
FeLoBALs may be associated with mergers. Further evi- 
dence comes from the discovery of iron absorption lines in 
the UV spectra of two low-redshift ULIRGs which harbor 
obscured AGN (iFarrah et al.ll2005D. Final ly, mid-IR pho- 
tometry observations (Farrah ct al. 2007a, hereafter F07) 
found that many FeLoBAL QSOs were IR-luminous, and 
may harbor high rates of star formation. 

Excellent distinguishing evidence of the 'youth vs. ori- 
entation' debate for FeLoBAL QSOs would be observa- 
tions that show directly that high rates of star forma- 
tion accompany the AGN. Given their reddened nature, a 
good place to look for such evidence is in the md-infrared, 
in which the spectral shapes of reddened AGN differ 
markedly from those of star-forming regions. The out- 
standing c apabilities of the Infrared Spectrograph (IRS, 
iHouck eta l. 2004) on-board the Spitzer space telescope 
(jWerner et al.l 2004) provided a huge step forward in 
available mid-IR spectroscopic capabilities, with the op- 
portunity to shed light on the FeLoBAL phenomenon. 
In this paper, we use Spitzer-IRS observations of six 
FeLoBAL QSOs to examine the idea that they are young 
objects. We assume a spatially flat cosmology, with 
Ho = 70 km s-i Mpc-\ 17 = 1, and ^a = 0.7. 

2. SAMPLE SELECTION 

We selected our sample in early 2006, with the re- 
quirements that candidates be confirmed as FeLoBAL 
QSOs via rest-frame UV spectroscopy, and lie in a red- 
shift range where we can observe useful diagnostics with 
the IRS. Accordingly, we chose six objects (Table [T]) at 
random from F07. The F07 sample is drawn randomly 
from the FeLoBAL QSO population known at the time, 
lie in the redshift range 1.0 < z < 1.8, thus placing im- 
portant mid-IR spectral features in the IRS bandpass, 
and have photometry observations at 24/ini, 70/ini and 
160/im, which provide good constraints on IR luminosi- 
ties. 

A downside of this selection though is that our sample 
is heterogeneous. When we selected our sample, the few 
known FeLoBAL QSOs had been found in several dif- 
ferent surveys. Since then, larger, homoge neous samples 
of FeLoBAL QSOs have been published (jTrump et al.l 
120061 ; iScaringi et al.l I2009D , but these samples were not 
available to us. Accordingly, some level of bias is in- 
evitable in our sample, though it is difhcult to quan- 
tify what effect this may have. We therefore sim- 
ply list the origins of our sample; th ree objects were 
found via the Sloan Digital Sky Survey (|Hall et al.ll2002l : 
lAdelman-McGarthv et al] [20081) . one (ISO 0056-2738) 
was disc overed serendipi tously from foUowup of distant 
clusters (|Duc et al.ll200l . and two (SDSS 1427-1-2709 & 
SDSS 1556-1-3517) were discovered during spectroscopic 



foUowup of quasars from the FIRST survey (jBecker et al.l 
[19971 INaiita"eraII[2000l) . SDSS 1556-H3517 is radio-loud, 
while the other five are radio quiet. 

Furthermore, at least three of our sam ple appear 
to ha ve atypicaO iron absorption features (jHall et al.l 
[200l . SDSS 1154-1-0300 has troughs where the ab- 
sorption remains significant at velocities comparable to 
the spacing between absorption features (so >12,000 km 
s^^), causing them to overlap each other. SDSS 2215- 
0045 and possibly SDSS 1214-0001 have stronger Felll 
absorption than Fell absorption, suggesting that the gas 
in which the BALs occur is dense, hot, and moderately 
highly ionized. We therefore do not correlate mid-IR 
spectral properties with those of the iron absorption fea- 
tures. 

3. METHODS 

We observed with the IRS using the first order of the 
short-low module, and both orders of the long-low mod- 
ule, giving observed-frame wavelength coverage of 7.5/xm 
to 35/^m. As the 24/^m fiuxes of our sample span a range 
of values, we used different observation times for each 
object so as to to give a signal-to-noise of at least 15 
in the continuum at observed-frame 24/^m. Observa- 
tions were performed in staring mode, using 'high' ac- 
curacy peak up observations performed with the blue ar- 
ray from a nearby Two Micron All-Sky Survey (2MASS, 
Ikrutskic ct al. 2006) star. 

The data were processed through the Spitzer Science 
Center's pipeline software (version 18.7), which performs 
standard tasks such as ramp fitting and dark current 
subtraction, and produces Basic Calibrated Data (BCD) 
frames. Starting with these frames, we produced re- 
duced spectra using the SMART v8.0 software package, 
follow ing the methods described in iLebouteiller et all 
()2010( ). which we summarize here. Individual frames 
were cleaned of 'bad' pixels using the IRSCLEAN task. 
The first and last five pixels, corresponding to regions of 
reduced sensitivity on the detector, were then removed. 
The individual frames at each nod position were then me- 
dian combined with equal weighting on each resolution 
element. Sky background was removed from each image 
by subtracting the image for the same order taken with 
the other nod position (i.e. 'nod- nod' sky subtraction). 
One-dimensional spectra were then extracted using 'op- 
timal' extraction with default parameters, and defringed 
using the internal SMART algorithm. We found that, 
in all cases, the sources were point-like, with a FWHM 
that was never wider than the PSF. We also checked 
that the sources were centered in the IRS slit by ex- 
tracting spectra using 'simple' (i.e. not PSF weighted) 
extraction. We found that the resulting spectra were of 
slightly lower signal-to-noise but consistent with the op- 
timally extracted spectra, confirming that any slit offset 
is insignificant. 

This procedure results in separate spectra for each nod 
and for each order. The spectra for each nod were in- 
spected; features present in only one nod were treated as 
artifacts and removed. The two nod positions were then 
combined. As the two nod positions have slightly offset 
wavelength grids, we combined the nods by first inter- 

^^ We include the qualifier as the rest-frame UV absorption line 
properties of FeLoBAL QSOs have not been studied exhaustively 
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leaving them, and then interpolating the fully sampled 
spectrum onto a reference wavelength grid. The nod- 
combined spectra in the three orders were then merged 
to give the final spectrum for each object. Overall, we 
obtained excellent continuum matches between different 
orders. Only in one case (SDSS 1556+3517) was there 
a significant order mismatch, between the SLl and LL2 
orders. For this object, we scaled the SLl spectrum by 
a factor of 1.10 to match the blue end of the LL2 contin- 
uum 

There remain however several uncertainties over er- 
ror propagation in the IRS reduction process. For ex- 
ample, defringing is still not completely understood, so 
some residuals likely remain that are not included in 
the 'formal' errors. Another example is the order mis- 
match between SLl and LL2 for SDSS 1556-f3517 - a 
PSF weighted 'optimal' extraction should in principle not 
produce this effect, and while the scaling is small, we do 
not understand why it is necessary. Overall therefore, we 
regard the resulting errors on each resolution element to 
be somewhat smaller than they should be, though the 
degree of this underestimate is likely insignificant. 

4. RESULTS 

The IRS spectra are presented in Figure [1] Spectral 
measurements are presented in Table [5] We found one 
further object serendipitously in the slit of ISO 0056- 
2738, which appears to be a PAH dominated system at 
z ~ 1.42, but do not consider it further in this paper. 

4.1. Spectral Features 

The spectra show a variety of spectral features. Four 
objects show one or more broad emission features at 
6.2/im, 7.7/im, 11.2/im and 12.7/im, attributed to bend- 
ing and stretching modes in Polycyclic Aromatic Hydro- 
carbons (PAHs, the 12.7/im feature also contains a con- 
tribution from [NeII]A12.81). The redshifts determined 
from the PAHs were in all cases consistent with the op- 
tical emission line redshifts, rather than the redshifts 
at which the UV absorption features peak, placing the 
source of the PAH emission in the host galaxy rather than 
in the outflow. The PAH fluxes and equivalent widths 
(EWs) were computed by integrating the flux after sub- 
tracting off a spline interpolated local continuum. An 
example of the spline fits is shown in Figure [51 

At least three objects contain a broad feature centered 
approximately at 9.7/xm, seen in both emission and ab- 
sorption, tha t arises from an Si-Q stre tching mode in 
Silicate dust (jKnacke fc ThomsonI 119731 ). We measured 
the strengths of these features via: 



. Fobs{9-7fim) \ 
^Fcont{Q-7f^m)J 



(1) 



where Fobs is the observed flux density at rest-frame 
9.7/im, and Fcont is the flux at the same wavelength de- 
duced from a splin e fit to the continuuni on either side 
(ISpoon et all 120071: ILevenson et all 120071: ISirockv et all 
l2008i r 

We also searched for other features seen in IR-luminous 
systems, though these are reliably measurable only with 
higher resolution data, so we do not present flux mea- 
surements. Both SDSS 1214-0001 and SDSS 1427-1-2709 
show strong [NeIII]A15.56, and weak but significant 



[NeV]A14.32. SDSS 1214-0001 additionally has weak de- 
tections at the positions of [Aril] A6.99 and H2S(2)A12.29. 
The other four objects show no further features that we 
can identify. SDSS 1154+0300 has an apparent feature 
at rest-frame ~ 8/im that we cannot find a reliable ID 
for, though as the spectrum is low S/N it is possible this 
feature arises from the juxtaposition of a declining PAH 
and a rising silicate feature. ISO 0056-2738 is too low 
S/N to infer the presence or otherwise of other features. 
We also note that all six spectra likely contain residual 
structure due to fringing that the defringing algorithm 
was unable to completely remove. 

4.2. Qualitative Comparisons 
4.2.1. Objects 3 & 4 

We start by comparing the two objects with promi- 
nent PAHs, SDSS 1214-0001 and SDSS 1427+2709, to 
well studied low-redshift objects (Figure [3|) . Both ob- 
jects closely resemble PAH dominated ULIRGs. It is 
superficially interesting that SDSS 1214-0001 is a good 
match to IRAS 15206+3342, a local ULIRG with iron 
absorption features in its UV spectrum, though not of 
great significance given that many local ULIRGs have 
similar mid-IR spectral shapes, see e.g. t he position of 
IRAS 15206+3342 in the 'network' plot of iFarrah et~all 
f2009b). SDSS 1427+2709 on the other hand is well 
matched by ULIRGs with weaker PAHs and a more pro- 
nounced continuuni, such as Mrk 231. Neither system 
resembles heavily absorbed ULIRG spectra such as Arp 
220. 

Moving on to comparisons with larger samples, we are 
hampered as there does not exist a comprehensive mid- 
IR spectroscopic survey of HiBALs or LoBALs to com- 
pare to. Even comparing to the general AGN population 
though, the peculiarity of these two objects, in partic- 
ular SDSS 1214-0001, is thrown into sharp relief. Their 
PAHs are extraordinarily prominent, far more so than for 
any Q SO, radio loud or radio quiet, that we are aware 
of (Haas et al. "20^ IHao et all [200l IShi et al.l [20061: 
iMaiolino et al. 200'<T), including those selected to be far- 
IR lu minous (jLutz et al.l 120081 : [Martinez- Sansigre et al.l 
l2008f). as well as the X-ray luminous 'type 2 QSO' objects 
(jSturm et al.l [20061. There are however a few systems 
with prominent PAHs among the Narrow Line Seyfert 
1 (NLSl, obje cts wi th H/3 FWHMs of < 2000km s-\ 
lOsterbrock fc Poggelll985D population (|Sani et al.l l20101. 
and several examples of Sy2 ULIRGs with strong PAHs 
(lArmus et al.l 120071 : llmanishi et al.l l2007t iFarrah et al] 
200^ see also the type 2 object LH901A in lSturm et al.l 
2d06h. 



4.2.2. Objects 1, 2, 5 & 6 

Turning to the other four sources; SDSS 2215-0045 
is a good match to QSOs with strong silicate emis- 
sion and weak but detec table PAHs, of which sev- 
eral exist (jHaoet al.l [20050 . Interestingly, SDSS 2215- 
0045 is an excellent match to PG1351+640, a IR- 
luminous, CO detected QSO with narrow BALs, Fell 
in emission un der H/3, and slight morphological distur- 
bance (Gclder nm^fc Whittle! Il99l iFalcke et al.l [19951 : 
IZheng et al...200lHEvans et al.ll200lD . Conversely, SDSS 
1556+3517, with its weaker silicate emission feature and 
negligible PAHs, is a good match both to some classi- 
cal QSOs, and to NLSl systems such as PG1211+143 
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(which has a hi gh velocity outflow, iPounds fc Pagel[2006l : 
iShi et alJ [2007h and IZwl (which has o ptical and UV 
Fell emiss ion and a decaying starburst (iMarziani et alJ 
[1996; Schi nnerer etHI [19981 : IS^race et al.lll998ll V SPSS 
1154+0300 has a poorer quality IRS spectrum but is also 
consistent with QSOs with weak silicate emission and 
negligible PAHs. ISO 0056-2738 is of too low signal-to- 
noise to draw any conclusions other than it appears to 
be consistent with AGN generally. 

4.2.3. The sample as a whole 

Finally, we compare the ensemble properties of our 
sample to other classes of active galaxies. We start 
by looking for matches with IR-selected AGN samples. 
One sample where we might e xpect ov erlap is with the 
X-ray detected LIRGs in .Brand et al.l ([2008b), but the 
tw o samples do not re semble each other. None of the 
16 iBrand et afl l2008bl objects have detectable PAHs; 
though 9 have some silicate absorption, the rest are fea- 
tureless power laws. Our two PAH dominated sources 
resemble some sources in op t ically faint TO/zm se lected 
samples ([Brand et al.l l2008al: [Farrah et al.l l2009a[ ). but 
70/zm selected samples have no sources with silicates 
in emission. We do better if we instead compare to 
IRS observations of high-redsh ift 15/xm selected samples 
([Hernan-Caballero et al.|[2009D . likely because we're se- 
lecting on hotter dust; Ibfim samples contain objects 
comparable to our PAH strong objects, but have few 
sources with silicates in emission. The one optically se- 
lected sample that ap pears reasonably matched to ours 
is the NLSl sample of iSani et al.l ([20101 ). which contains 
both PAH dominated objects, and a few objects with 
what appears to be weak silicate emission. 

4.3. Spectral Diagnostics 

We move on to quantitative diagnostics. As our spec- 
tra are of relatively low S/N and do not have coverage 
in all the IRS modules, we use simple diagnostics that 
allow for easy comparisons with other samples. 

We start with the PAH features. The PAH flux ratios 
of our sample, considered either as functions of each other 
or of PAH luminosity (Figures [4| and [5|) are slightly offset 
from those of low-redshift ULIRGs, but lie within the dis- 
persion of high-redshift 70/im or 24/r selected samples. 
This is straighforward to understand. Local ULIRGs are 
selected without a bias towards AGN, and have lower IR 
lum inosities, on average, t han our samp l e. Co nversely, 
the iSaiina et al] (|2007[i & iDasvra et all ([2OO90 samples 
are (arguably) biased towards AGN, and have compara- 
ble total IR luminosities to our objects. It is thus not sur- 
prising that the PAH flux ratios and luminosities of our 
sample resemble the high redshift comparison objects. It 
is also interesting that the two systems with the strongest 
PAH detections (SDSS 1214-0001 & SDSS 1427-1-2709) 
are also the strongest 160fim detections, consistent with 
the idea that starbursts are associated with colder dust, 
but the small size of our sample means this consistency 
could simply be coincidence, so we do not comment on 
it further. 

We estimate star formation rates from the PAH fea- 
tures using: 



where Lp is the combined l uminosity of the 6.2/xm 
and 11.2/.tni PAH features (Fa rrah et al]|2007b[ ). The 
values are listed in Table [2PI . Overall, they are 
comparable to those derived for other high redshift 
IR-lumi nous sources , incl uding 'bump' selected star- 
bursts ([Farrah et "all 120081), sub-mm selected starbusts 
([Pope et al.ll2008HMenendez-Delmes t re et al.l[200a\ and 
far-IR bright OSOs (|Lutz et al.l 120081 ). The star forma- 
tion rate in SDSS 1214-0001 is extraordinarily high; as- 
suming continuous star formation then this star forma- 
tion rate is capable of manufacturing IO^^Mq of stars in 
less than lOOMyr, and < 50Myr if we assume a late stage 
exponentially decaying burst. 

Ne xt, we employ the 'Fork' diagnostic of iSpoon et al.l 
120071 (Figure[6]), which employs both the 6.2^m PAH fea- 
ture and the 9.7 fim silicate feature. Here our sample is 
not distributed in a similar way to the 70/im or 24/r se- 
lected so urces, but instead l ies along the lower branch of 
the fork. ISpoon et al.l 120071 postulate that sources move 
around the Fork diagram as their power source evolves, 
starting on the upper branch and then moving diago- 
nally or vertically downwards as starburst/AGN activity 
clears the obscuration from the nuclear regions. There- 
fore, solely from this diagnostic, we would classify our 
sample as late-stage ULIRGs. 

We move on to consider mid-IR continuum diagnos- 
tics. Given the limited wavelength coverage, we use 
the rest-frame 6/im continuum lur ninosity a s a proxy 
for the bolometric A GN luminos itv (iNardini e t al. 200l; 
iWatabeet aI[|2009D . Following INardini et al.l.2008. we 
define: 



R 






(3) 



where Lq is the 6/im luminosity as measured from the 
IRS spectra and Lm is the 1-1000/im luminosity from 
F07. The fractional contribution of an AGN to the 6/im 
luminosity, ae, is then: 



1 



ae 



f RaRs 



Rs~Ra\ R 
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(4) 



where Rs and Ra are the equivalents of R for 'pure' 
starbursts and AGN, wi th values of (117 ib 8 ) x lO"'^ and 
0.32 ±0.1 respectively (INardini et al.ll2008[) . Therefore, 
the fractional contribution of the AGN to the total IR 
luminosity, Uboi, is: 



aboi 



ae 



ae 



(l-ae) 



(5) 



SFR[MQyr- 



1.18 X IQ-'^'^ Lp[ergs s"^] (2) 



The resulting atoi values are listed in Table [T] For com- 
parison, we also list the atoi values obtained from the 
SEDs in F07. In four cases the values are consistent, 
but for two (SDSS 1556-^3517 & SDSS 2215-0045) they 
are not. We think it likely that the SED derived val- 
ues are more reliable, as they are direct measurements 
from the (albeit sparsely sampled) full SEDs, while the 
values computed using Equation [S] are calibrated using 
a wide range of sources. It is interesting that these two 

^^ Estimating t he star formation rates using the formula in 
IHouck et aLI l|2007l ) gives comparable results 
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objects are the only two with detected sihcate emission, 
but we do not know if this is the cause of the discrep- 
ancy, or coincidence. Considering either method though, 
the sources with strong PAHs have weak AGN contri- 
butions to the total IR luminosity, and the range in aboi 
values for the whole sample is similar to that seen in local 
ULIRGs. 

Finally, we fit the source with the mos t prominent 
PAH s, SDSS 1214-0001, with PAHFIT (|Smith et alj 
|2007| ). This is not an attempt to reconstruct star for- 
mation parameters as PAHFIT is intended for systems 
where emission from stars provides at least most of the 
flux across the mid-IR, but it does serve as a test of 
how much of the mid-IR flux in SDSS 1214-0001 is at- 
tributable to star formation. The result is shown in Fig- 
ure [T] The fit is good, with xled — l-^- The fit is poor at 
< 6/im, with a steeply rising contribution from 'starlight' 
with decreasing wavelength that is probably an attempt 
to fit the AGN continuum. At A > 6/im though, we ob- 
tain an excellent fit. The PAHs are well fitted by the 
Drude profiles, and the continuum is predicted to mostly 
come from dust at 50k, with small contributions from 
the 135k and 300k components. 

5. DISCUSSION & CONCLUSIONS 

A general caveat to all of what follows is our small and 
heterogeneous sample. Our conclusions should thus be 
regarded as tentative. 

The mid-IR spectra of our sample span a wide range 
of shapes. We see classical QSO spectra with hot sili- 
cate dust, together with classical starburst spectra with 
strong PAHs. The 'starburst' spectra have more promi- 
nent PAHs than those seen in any other QSO so far ob- 
served, while the hot silicate dust sources do not have 
a close match in any purely mid/far-IR selected sample 
that we are aware of. It is possible that our heteroge- 
neous selection leads to the heterogeneity of the spectra, 
but if this were true then we might expect the ISO se- 
lected object to have the strongest PAHs, which is not 
the case. Indeed, the two spectra at the extreme ends 
of the range of spectral shapes are both SDSS objects. 
The star formation rates of our sample span levels com- 
parable to those in the most luminous starbursts found 
in any survey at any wavelength, to those seen in moder- 
ately IR-luminous AGN. The spectral diagnostics mark 
our sample as late-stage ULIRGs with a wide range of 
fractional AGN contributions. 

These results, combined with their optical classifica- 
tion as reddened QSOs with a strong outflow, are in prin- 
ciple compatible with a continuum of 'end-of-ULIRG' vs 
'peculiar QSO' contributions to the FeLoBAL QSO pop- 
ulation. To frame the discussion, we describe three points 
on this continuum in detail: 

1: FeLoBAL QSOs are a transition stage between star- 
burst dominated ULIRGs and QSOs. As our sample 
span the entirety of this transition, and because FeLoB- 
ALs are observed to be rar43i we can set constraints 
on the timescale of this transition depending on how in- 
trinsically common the Fe absorption features are. If Fe 

^^ A further reason behind the rarity of FeLoBALs is that they 
they are hard to detect; FeLoBAL QSOs are red, have few or no 
broad emission hnes and no UV ex cess, and are thus hard to find 
in standard QSO searches (see e.g. lAppenzeller et 311120051 ). The 
iron absorption features themselves are weak. 



absorption is ubiquitous in such a transition, then the 
transition must be short in comparison to the ULIRG or 
QSO hfetime. Ta king the ULIRG a nd QSO lifetimes to 
be '-^lO* years (T acconi et al . 1 120081 ) . then the transition 
must take ^10^ years (see also lGregg et al.ll2002[ ). If on 
the other hand the iron absorption features are intrinsi- 
cally rare, then the transition can take longer. 
2: FeLoBALs are comprised of comparable fractions of 
two galaxy populations, one that is transitioning from a 
ULIRG to a QSO, and another that is an unusual phase 
that QSOs go through, unconnected to mergers. Here 
the heterogeneity of our sample arises from observing 
two unconnected populations. 

3: FeLoBAL QSOs are entirely an 'unusual QSO' class, 
and have nothing to do with mergers or otherwise 'tran- 
sitioning' systems. Instead, the iron absorption features 
mark QSOs with both an outflow and atypically large 
amounts of iron in their ISM. 

Formally, we cannot rule any of these scenarios out. 
The third though seems unlikely. It is hard to see how 
such a heterogeneous set of mid-IR spectra could arise 
from observing a peculiar class of QSO in which nothing 
more interesting than an outflow with large quantities 
of iron is going on, especially if we assume that the re- 
strictive selection on rest-frame UV properties means we 
are viewing them over a restricted range in line of sight. 
Furthermore, most systems with star formation rates in 
the range seen in our sample are ULIRGs, which (at low 
redshift, at least) are nearly all mergers. 

Distinguishing between the first and second scenarios 
is however more difficult. On one hand, the first scenario 
uses a single origin for which there is independent evi- 
dence from studying ULIRGs. The second scenario thus 
seems contrived, as it proposes two origins where one 
will do. On the other hand, our sample is small, and it 
is tempting to place undue emphasis on the properties of 
SDSS 1214-0001, as it is so peculiar. If this object were 
removed then the 'unusual QSO' scenario would become 
more attractive. That said, the counter-argument also 
holds; remove (say) SDSS 1556-1-3517, and the 'end-of- 
ULIRG' scenario becomes more attractive. We are mind- 
ful though that it is only the properties of SDSS 1214- 
0001 (and to an extent SDSS 1427-f 2709) that make us 
seriously consider the 'end-of-ULIRG' scenario. 

Recent work that may s hed li ght on this is the study 
of NLSl's by iSani et al.l ([2010). Their spectra show 
that at least some NLSl's harbor intense star forma- 
tion and a weak mid-IR AGN continuum. NLSl's also 
have strong optical Fe emission lines, a large soft X-ray 
excess, and exhibit rapid, large amplitude X-ray vari- 
ability. Furthermore, their black hole masses appear to 
be smaller than those in BLSl's of comparable luminos- 
ity, tho ugh there is controversy over whether NLSl's lie 
bel ow (iPeterson et all l2000t iGrupe fc Mathuii I2004D or 



on (iBotte et al.ll2005l : lKomossa fc Xull2007l:lDecarn et al.l 
|200^ the7lf^-cr* relation (jTremaine et al.ll2002D . On 
one level this hints at interesting links between NLSls 
and FeLoBAL QSOs; for example, FeLoBAL QSOs may 
be analogues of NLSls seen more pole-on than edge-on, 
and in which a strong outflow has formed. Exploring 
this idea in detail requires significantly larger samples, 
so we do not pursue it here. More generally however, it 
implies that relative orientation may also play a role in 
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the FeLoBAL QSO phenomenon. 

We think it hkely that the simplest solution that is 
consistent with all the observations, and the apparent 
rarity of FcLoBALs, is the correct one. We therefore pro- 
pose that FeLoBAL QSOs as a class are rapidly evolving, 
youthful QSOs in which a starburst is coming to an end, 
and the AGN is in the last stages of burning through its 
surrounding dust. We also propose, with more reserve, 
that (1) we view FeLoBAL QSOs over a restricted line 
of sight range, where we see them more pole-on than 
edge on, and (2) the outflow has in some cases partially 
cleared the dust from around the AGN, leaving the IR 
emission from the starburst to dominate. 

There are three ways to test this conclusion. First, 
FeLoBAL QSO hosts should be in the final stages of 
merging, so their host galaxies should show slight signs of 
morphological disturbance. Second, larger mid-IR spec- 
troscopic surveys of FeLoBAL QSOs should show a simi- 
lar range of spectral shapes to ours. Third, far-IR photo- 
metric surveys of FeLoB ALs should find both a moderate 
enhancement of the fraction of far-IR bright sources com- 
pared to the classical QSO population, corresponding to 
those FeLoBAL QSOs that still harbor high star forma- 
tion rates, and a correlation between far-IR luminosity 
and PAH equivalent width. Moreover, it would be useful 
to perform radiative transfer modelling of a BAL wind 
through a starburst to determine the conditions under 
which iron absorption is observed in such systems, and 
so provide insight into whether our conclusions are fea- 
sible. 

We close with a caveat and a speculative comment. 
First, there is an important parameter that we cannot 
address; the way in which iron absorption features in 
FeLoBALs signpost starburst and AGN activity. For ex- 
ample, if the Fe absorption only occurs when a starburst 
is ending and an AGN is blowing away its surround- 
ing dust, but can occur throughout such a transition, 
then our first scenario is likely correct. As a contrived 
counter-example, if Fe absorption can only occur dur- 
ing the early phases of such a transition, but can occur 
randomly in reddened AGN in which there is no signifi- 
cant star formation, then in a small sample such as ours 
we would conclude, incorrectly, that all FeLoBAL QSOs 



were a ULIRG-to-QSO transition where the Fe absorp- 
tion is present throughout. Resolving this issue is beyond 
the scope of this paper, so we note it as a caveat to our 
conclusio ns. ^_^ 

Second, Hall et aT (2002) note that two of our sample, 
SDSS 1214-0001 and SDSS 2215-0045 have [Felll] / [Fell] 
ratios greater than unity, implying the BALs arise in 
unusually dense, hot gas. The two ratios however 
are different; SDSS 2215-0045 has a [Felll] / [Fell] ratio 
much greater than unity, while SDSS 1214-0001 has an 
[Felll] / [Fell] ratio only slightly greater than unity. While 
conclusions drawn on only two objects are not trustwor- 
thy, this difference in ratio is consistent with our pro- 
posed evolutionary sequence. SDSS 1214-0001 has strong 
PAHs and sihcates in absorption, while SDSS 2215-0045 
has weak PAHs and silicates in emission, implying that 
SDSS 1214-0001 is at an earlier evolutionary phase than 
SDSS 2215-0045. Moreover, the enhanced [Felll] / [Fell] 
ratio in SDSS 2215-0045 compared to SDSS 1214-0001 
could be interpreted as the BAL outflow in SDSS 2215- 
0045 being more developecP^. Further support for this 
idea would be an [Felll] / [Fell] ratio greater than unity 
in SDSS 1427-1-2709, but the only publi shed spectrum 
we are aware of is in I Becker et al.l (|2000( ) , in which the 
Felll lines lie in a noisy region at the extreme blue end 
of the bandpass. They do not however appear to be sig- 
niflcantly stronger than the Fell lines, so we simply note 
this as an interesting avenue for future work. 
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TABLE 1 
FeLoBAL QSO Sample 



Object 



RA (2000) 



Dec 



^sys 



^abs 



/to /i60 Ofcoilefjm Ctbolls 



(1) ISO J005645. 1-273816 

(2) SDSS J115436. 60+030006.3 

(3) SDSS J121441. 42-000137.8 

(4) SDSS J142703.62-f 270940.3 

(5) SDSS J155633. 78+351757.3 

(6) SDSS J221511. 93-004549. 9 



00 56 45.2 


-27 38 15.6 


1.78 


1.75 


1.6 


<7.0 


<50 


0.09- 


■ 0.40 


>0.26 


12.5-13.0 


11 54 36.6 


03 00 06.4 


1.46'= 


1.36 


7.6 


17.9 


<50 


0.35- 


■ 0.60 


>0.40 


12.9-13.1 


12 14 41.4 


-00 01 37.9 


1.05 


0.99 


4.7 


38.3 


78.9 


0.07- 


■ 0.15 


<0.37 


12.7-12.9 


14 27 03.6 


27 09 40.3 


1.17 


? 


4.8 


32.6 


68.1 


0.09- 


■ 0.17 


<0.61 


12.8-13.0 


15 56 33.8 


35 17 58.0 


1.50 


1.48 


13.9 


24.7 


<50 


0.40- 


■ 0.70 


>0.75 


13.1-13.3 


22 15 11.9 


-00 45 49.9 


1.48 


1.36 


10.4 


27.2 


<50 


0.09- 


■ 0.30 


>0.36 


13.0-13.4 



Note. — All flux densities are quoted in mjy. Errors are typically 10% at 24/^m, 20% at 70fini, and 25% at 160/xni. Luminosities are the 
logarithm of the rest-frame 1-1000/im luminosity, in units of solar luminosities (3.826x10 Watts), taken from F07. Limits and ranges are 

3CT. 

^ Systemic rcdshift from narrow optical emission linos 

Redshift of peak absorption from the broad UV absorption lines 
^ IHall et J1 (i200^ 

Fractional AGN luminosity, computed using the prescription in ^4.31 ( iNardini et al.||2008f) . 
^ Fractional AGN luminosity, derived from the SED fits in F07. 



TABLE 2 

Spectral Measurements 



Object 


PAH 6.2/^m 


PAH 7.7/im 


PAH 11.2^111 


Ssil 


SFR^' 




Flux 


EW 


Flux 


EW 


Flux 


EW 




Mq yr-i 


(1) ISO 0056-2738 


6.1+8.9 


0.06+0.09 


6.0+8.0 


0.08+0.12 


4.1+3.9 


0.08+0.08 


0.25+0.11 


2600 + 2500 


(2) SDSS 1154+0300 


5.2+5.7 


0.01+0.01 


15.1+9.2 


0.04+0.02 


0.3+4.6 


0.01+0.02 


0.15+0.08 


900 + 1100 


(3) SDSS 1214-0001 


16.3+3.3 


0.06+0.01 


56.5+7.0 


0.25+0.05 


22.5+5.6 


0.24+0.06 


-0.25+0.05 


2700 + 500 


(4) SDSS 1427+2709 


7.3+2.6 


0.03+0.01 


25.0+6.4 


0.10+0.03 


2.4+2.1 


0.02+0.01 


-0.10+0.05 


900 + 350 


(5) SDSS 1556+3517 


1.0+3.0 


0.01+0.01 


11.2+6.6 


0.02+0.01 


3.0+3.9 


0.01+0.02 


0.22+0.04 


700 + 800 


(6) SDSS 2215-0045 


10.5+5.3 


0.03+0.01 


28.1+6.8 


0.09+0.02 


7.7+3.7 


0.02+0.01 


0.62+0.10 


2700 + 1000 



Note. — PAH fluxes arc in units of 10 W cm and equivalent widths in fim. 

Star formation rate, determined from Equation[2] Error is solely that derived from the uncertainty in the fluxes. 
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Fig. 1. — Spitzer-IRS low-resolution spectra of our sample, plotted in the rest-frame using the optical emission line redshifts. The data 
are plotted in black, with \u errors in grey. Flux densities are in injy and wavelengths are in /^m. The number in the top left of each panel 
is the ID number in Table[T] The green solid (dashed) lines mark the 6.2/im, 7.7/^m, 8.6/^m, W.l^va and 12.7/.tm PAH features assuming 
the systemic (peak absorption, if available) redshift. The red lines perform the same function for the [Ne IIJA12.81, [Ne V]A1 4.32 and [Ne 
III]A15.56 li nes. The yellow shaded region shows the approximate extent of the 9.7/im silicate absorption feature, but see IBowev et all 
[T998. : ,Draini[200a : INiicutta et al.„200a . 



10 



Farrah et al 




5.5 



6.0 6.5 

Rest wavelength (urn) 



7.0 



Fig. 2. — An example of the spline fits used to determine the shape of the underlying continuum for the PAH fluxes, in this case for the 
6.2/^m PAH feature in SDSS 1214-0001. The (rest-frame) spectrum is plotted in black. The points used to define the continuum are shown 
in green, while the spline fit and its residuals are shown in blue. The errors on the spectrum have been omitted for clarity. The dashed line 
shows the adopted peak wavelength of the PAH feature, while the dotted lines show the lower and upper limits of the integration. This 
plot also demonstrates the excellent agreement between the two nods; the 'double' blue error bars show separately the error bars for the 
two nod positions, which have slightly offset wavelength grids. 
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Fig. 3. — Comparison of our spectra to well studied low-redshift obj ects ^^_^ ^_^ .».^^^_^^^_ 

the others are in the right panel. Comparison spectra are taken from iBrandl et al.l l|2006l ) : lArmus et al.l 
llma nishi et al. (2007). Left: Red - Mrk 231. Orange - composite of local ~ IO^^Lq starbursts. Green - Arp 220 
(a local ULIRG with iron absorption features in its UV spectrum). Right: Red - 3C 273. Orange - I Zw 1 (a low 
1211-1-143. Blue - PC 13514-640 (both far-IR luminous QSOs). 



t panel, and 

I); IShi et al.l pOOTT l and 

Blue - IRAS 15206-1-3342 

NLSl). Green - PG 
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Fig. 4. — PAH ratio diagnostic plots, divided by comparisons t o (top row ) low-rcdshift and (bottom row) high redshift samples. Black: 
FeLoBALs. Green: low-z ULIRGs with measurable 7.7/im PAHs l|SDOon et al. 2007; Dcsai et al. 2003)- Red: low-z AGN (Wccdman ct aQ 



09al ) . Grey: high-z 'bump' selected (Farrah ct al. 2008). Orange: high-24/im/0.7^m selected Dasvra et al. (2009). Cyan: high-z 
24um/8/im and 24/im/0.7/im selected iSaiina et al.i (i2007i ). Yellow: high-z sub-mm selected {Pope et al.ii2008i : iMenendez-Delmestre et al.l 
[2009D . 
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Fig. 5. — PAH luminosity diagnostic plots. Color coding is the same as in Figure HI 
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Fig. 6. — 'Fork' diagnostic diagram dSpo on ct al."2007). The FeLoBALs arc plotted in black, while the other points are color-coded as 
in Figure |4] Plots of the composite spectra for each region are in figure 2 of Spoon et al. (2007). With respect to the comparison objects 
plotted in Figure [S] Mrk 231 is class lA, IRAS 15206+3342 is class IB, Arp 220 is class 3B, the starburst composite is class IC, and the 
remaining objects are all class lA. The 'bump' sources lie close together on this plot, so, for clarity, we plot their average as a single point. 
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Fig. 7. — PAHfit results for SDSS 1214-0001. Details of the model parameters can be found in lSmith et al.l II2007I '). Vertical black lines: 
observed data. Green: combined best-fit model. Red solid lines: thermal dust continuum components. Magenta: starlight continuum. 
Black solid: combined starlight+thermal dust continuum. Violet: fine structure+molecular features. Blue: PAH features. Black dotted: 
extinction curve. 



